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A Rational Approach to the Resonance Region 
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Resonance Saturation in QCD can be understood in the large- N c limit from the mathematical theory of Pade 
Approximants to meromorphic functions. 



in: 



> 
O 

(N 

On 
O 
OO 
O 



X 



1. Introduction 

The Chiral Lagrangian organizes the physics 
of the strong interactions at low energy as an 
expansion in powers of momentum and masses 
of the lightest pseudoscalar fields. Since all the 
heavier states of QCD are integrated out, their 
physics is encoded in a set of Low Energy Con- 
stants (LECs). These LECs are indispensable 
to make definite predictions in Chiral Perturba- 
tion Theory. However, while at 0(p 4 ) the LECs 
are relatively well known, at 0(p 6 ), most of the 
0(100) LECs are completely unknown. 

At the same time, QCD Green's functions seem 
to be approximately saturated by a few reso- 
nances. In the vector channel, that fact was called 
Vector Meson Dominance. Also, LECs seem to 
be well saturated by the lowest meson (after con- 
straints in the OPE) [2]. 

The theoretical support for these phenomeno- 
logical ideas comes out from the framework of 
Large- N c QCD. The main point is that in Large- 
N c QCD mesons are stable [3]. With this con- 
sideration and keeping only a finite set of reso- 
nances (instead of an infinite number of mesons 
in Green's functions), [1] proposed the Minimal 
Hadronic Approximation (MHA). The question 
now is that MHA only uses a finite set of reso- 
nances. How the error can be estimated and how 
the approximation can be improved is the main 
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goal of [T] within the mathematical theory of Pade 
Approximants. 

Let f(z) be a function of a complex z with 
and expansion f(z) = J2^=o fn,z n when z — > 0. 

'n 



Then, we can define a rational function P^f (z) 



such that PM(Z) EE « /„ + hz + f 2 z 2 
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Then, P, 



N ( z ) 



• ■ • + Jm+nz ^ + U(z ) 
is a Pade Approximant (PA). 

The convergence properties of the PAs to a 
given function are much more difficult than those 
of normal power series. However, those which 
concern meromorphic functions are rather well 
known. The main result which we will use is 
Pommerenke's Theorem [5] which asserts that the 
sequence of PAs to a meromorphic function is 
convergent everywhere in any compact set of the 
complex plane except, perhaps, in a set of zero 
area. This set includes the set of poles where the 
original function f(z) is ill-defined but there may 
include extraneous poles as well. The previous 
convergence theorem requires that, either these 
extraneous poles move very far away from the re- 
gion as the order of the PA increases, or they pair 
up with a close-by zero (Fig. [1]) becoming what 
is called a defect [6]. These must be considered 
artifacts of the approximation. Near these extra- 
neous poles the approximation breaks down, but 
away from them, the approximation is safe. 

In the physical case, the original function f(z) 
will be a Green's function of the momentum vari- 
able Q 2 . In Large- N c QCD this Green's function 
is meromorphic with all its poles located on the 
negative real axis in the complex Q 2 plane. Its 
poles are identified with the meson masses. On 
the other hand, the region to be approximated 
by the PAs will be that of euclidean values for 
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the momentum Q 2 > 0. The expansion of the 
Green's function for large and positive Q 2 coin- 
cides with the Operator Product Expansion. We 
will see that some of the poles and residues of the 
PAs may become complex. The reason lies in the 
fact that a meromorphic function does not obey 
any positivity constraints. This means that these 
poles and residues may have nothing to do with 
the physical meson masses and decay constants. 
However, and this is pointed out in pQ, this does 
not spoil the validity of the rational approxima- 
tion provided the poles, complex or not, are not 
in the region of Q 2 one is interested in. It is to be 
considered rather as the price to pay for using a 
rational function, which has only a finite number 
of poles, as an approximation to a meromorphic 
function with an infinite set of poles. 

2. The model 

To exemplify these characteristics, let consider 
the two-point functions of vector and axial- vector 
currents in the chiral limit 



= * J d*xe i «*(jV> A (x)4 v > A (0))= (1) 
{q^iu - g^q 2 )Rv,A(q 2 ) , 

with Jy(x) — d(x)j fl u(x) and J A {x) = 
d(x)j^j 5 u(x). Following Refs. [10], we define 
our model as {Tl LR (q 2 ) = ^(U v (q 2 ) - U A (q 2 ))): 
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q 2 -q 2 + M 2 
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-q 2 + M 2 (n) -q 2 + M 2 A (n) 



Here F p , M p are the electromagnetic decay con- 
stant and mass of the p meson and Fy^in) are 
the electromagnetic decay constants of the n — th 
resonance in the vector (resp. axial) channels, 
while Mv,A{n) are the corresponding masses. Fq 
is the pion decay constant in the chiral limit. The 
dependence on the resonance excitation number 
n is Fy A (n) = F 2 = constant and M v A (n) — 
my A + n A 2 in accord with known properties of 
the large- N c limit of QCD as well as alleged prop- 
erties of the associated Regge theory [5] . 



To resemble the case of QCD, we will demand 
that the usual parton-model logarithm is repro- 
duced in both vector and axial-vector channels 
and that TlL R (q 2 ) has an operator product expan- 
sion which starts at dimension six. A set of pa- 
rameters satisfying these conditions, while keep- 
ing the model realistic and at a manageable level, 
is given by F Q = 85.8 MeV, F p = 133.884 MeV, 
F = 143.758 MeV, M p = 0.767 GeV, m A = 
1.182 GeV, my = 1.49371 GeV and A = 
1.277 '4 GeV. The values are chosen so that: 
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2.1. PA with a Model 

The simplest PA to the function Q 2 T\ LR (-Q 2 ) 
with the right fall-off as Q~ A at large Q 2 is 
P^Q 2 ), which has three unknown coefficients 
fixed with the first three coefficients from the chi- 
ral expansion of ([3]), i.e., Co.2,4: 



P 2 °(Q 2 ) 



(Q 2 + z R ){Q 2 + z* R ) 



(5) 



where r% = 3.379 x 10~ 3 and z R = 0.6550 + 
i 0.1732. 

A pair of complex-conjugate poles means that 
these poles cannot be interpreted in any way as 
the meson states appearing in the physical spec- 
trum. In spite of this, if one expands ([5]) for large 
Q 2 > 0, one finds CU = -r% = 3.379 x 10" 3 
which differs from the real value by about 30%. 
Even better is the prediction of the fourth therm 
in the chiral expansion (3% of error). This agree- 
ment is not a numerical coincidence and the ap- 
proximation can be systematically improved if 
more terms of the chiral expansion are known 
(Fig. [TJtop). For example, going up to P^ 
(with 103 parameters), one finds that this ratio- 
nal approximant correctly determines the values 
for C_4._6,-8 with 52,48,45 decimal figures re- 
spectively. In the case of C2061 which is the first 
predictable term from the chiral expansion, the 
accuracy reaches 192 decimal figures. 

As it happens for ([5]), also higher order PAs 
may develop some artificial poles. Fig. [TJbottom 
shows the location of the 52 poles of the PA in 
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Figure 1. Top: Location of the poles of P$ (green), 
P32 (blue) and P52 (red) in the complex q 2 plane. 
When the order of the PA is increased, the over- 
all shape of the figure does not change but the two 
branches of complex poles move toward the right, 
away from the origin; Bottom: Location of the poles 
(dots) and zeros (squares) of P52 in the complex q 2 
plane. 



the complex q 2 plane. A detailed numerical anal- 
ysis reveals that the poles and residues reproduce 
very well the value of the meson masses and de- 
cay constants for the lowest part of the physical 
spectrum of the model given in ((21) , but the agree- 
ment deteriorates very quickly as one gets farther 
away from the origin, eventually becoming the 
complex numbers seen in Fig. [T] It is by creating 
these analytic defects that rational functions can 
effectively mimic with a finite number of poles 
the infinite tower of poles present in the original 
function ([2]). Then, the determination of decay 
constants and masses extracted as the residues 
and poles of a PA deteriorate very quickly as one 
moves away from the origin. There is no reason 
why the last poles and residues in the PA are to 
be anywhere near their physical counterparts and 



their identification with the particle's mass and 
decay constant should be considered unreliable. 
Clearly, this particularly affects low-order PAs. 

3. A QCD case 

In the case of large- QCD in the chiral limit, 
two considerations arise. First, the lack of input 
data makes going to higher order in the construc- 
tion of the rational approximants a difficult work. 
In addition, an order of PAs is increased by one 
unit by adding two new inputs. Second, any in- 
put value will have an error (both experimental 
and because of chiral and large- N c limits), and 
this error will propagate through the rational ap- 
proximant. In spite of these difficulties one may 
feel encouraged by the phenomenological fact that 
resonance saturation approximates meson physics 
rather well. 

The simplest PA to the function Q 2 U LR (-Q 2 ) 
with the right fall-off as Q~ 4 at large Q 2 is 

P^iQ 2 ) — i+A q°+-B q 4 • T nc vaiucs °f the three 
unknowns a, A, B may be fixed by requiring that 
his PA reproduces the correct values for F ,Li 
and SM^, the electromagnetic pion mass differ- 
ence (which is defined in terms of U LR as 5M 2 = 
M 2 + -M 2 = -^J^dQ 2 Q 2 U LR (Q 2 )). 
Then, we use: F = 0.086 ± 0.001 GeV, SM n = 
4.5936 ± 0.0005 MeV, L w {0.5 GeV) < L w < 
L w (l.l GeV) L w = (-5.13 ±0.6) x 10" 3 . 

The final results are shown in table [TJ In 
a real QCD case, the lack of input data does 
not allow us systematically improve our results. 
However, other Pade Approximants are useful for 
doing this improvement, for instance, what is 
called Pade- Type Approximants (PT). PTs are 
defined as the ratio T^f of two polynomials where 
the polynomial in the denominator has its zeros 
given in advance, preassigned, for example, at 
the position of the first N resonance masses on 
the original Green's function. Therefore, we will 
use for the masses [9]: m p = 0.7759 ± 0.0005, 
m p > = 1.459 ± 0.011, m p „ = 1.720 ± 0.020, 
m p n, = 1.880 ± 0.030,m ai = 1.230 ± 0.040 and 
= 1.647 ± 0.022 (all numbers in GeV). Be- 
cause we can compute succesive terms in a PT se- 
quence, we can get a feeling of convergence and es- 
timate in this way the error of our approximation. 
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-43 ± 13 

— 86*87 


81 ±53 


-145 ± 120 


-4.1 ±0.5 


6±2 


-7 ±6 



Table 1 



Values of the coefficients C2k in the high- and low-Q 2 expansions of Q 2 Hlr(—Q 2 ) in Eq. in units of 

1(T 3 GeV 2 ~ 2k . Recall that C- 2 = 0. 



Our final result is C 87 = (5.7±0.5) x IQ^GeV 2 . 
Recently, [7] has found C 87 (m = 0.77GeV) = 
(3.9 ± 1.4) x lO- 3 GeV~ 2 at NLO in l/N c in Res- 
onance Chiral Theory. 
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Figure 2. Prediction for C«t in the large- A^c limit 
from the PA P° an d PTs in Ref . |Tj . For comparison 
we also show the estimate from Refs. [8j, which we 
label 'A', 'B' and 'C (resp.). 



4. Conclusion 

Ref. [1] pointed out that approximating Large- 
N c QCD with a finite number of resonances may 
be reinterpreted within the mathematical Theory 
of Pade Approximants to meromorphic functions. 
With the help of a toy model, we have reviewed 
the main results of this theory: rational approxi- 
mants to a meromorphic Green's functions yield 
an accurate description of the Green's function 
in the Euclidean region Re(q 2 ) < while creating 
artificial poles (and corresponding residues) in the 
Minkowsi region Re(q 2 ) > 0. Therefore, althoug 
they may be a very good approximation in the eu- 
clidean, it is in general unreliable to extract prop- 
erties of individual mesons, such as masses and 



decay constants from an approximation to Large- 
N c QCD with only a finite number of states. 

Pade Approximants beyond Large-A^ c QCD 
have been applied, recently, to the study of the 
pion form factor pTj . 
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